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Properties of the Kondo effect in quantum dots depend sensitively on the coupling parameters and so on the 
realization of the quantum dot - the Kondo temperature itself becomes a mesoscopic quantity. Assuming chaotic 
dynamics in the dot, we use random matrix theory to calculate the distribution of both the Kondo temperature 
and the conductance in the Coulomb blockade regime. We study two experimentally relevant cases: leads with 
single channels and leads with many channels. In the single-channel case, the distribution of the conductance 
is very wide as Tk fluctuates on a logarithmic scale. As the number of channels increases, there is a slow 
crossover to a self-averaging regime. 
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Advances in fabrication of nanoscale devices have made 
possible unprecedented control over their properties. 1 While 
the transport properties of quantum dots - systems of electrons 
confined to small regions of space - have been studied exten- 
sively for the last decade, it is only recently that their many- 
body aspects have been probed 2,3 using the exquisite control 
now available. On the other hand, a continually fascinating 
aspect of nanophysics is the presence of quantum coherence 
and the interference "fluctuations" that it engenders.^ It is 
natural, then, to ask how this classic mesoscopic physics ef- 
fect influences the newly found many-body aspects. 

Quantum dots are usually formed in one of two ways: ei- 
ther by depleting a two-dimensional electron gas (2DEG) at 
the interface of a GaAs/AlGaAs heterostructure - a semicon- 
ductor quantum dot (SQD)I - or by attaching leads to ultra- 
small metallic grains (MQD). 2 If the dot is weakly coupled 
to all leads used to probe it, there is an energy price to be 
paid for an excess electron to enter the dot, simply because 
of the energy required to localize the charge. This regime is 
known as the Coulomb Blockade (CB). In a SQD, one attains 
this regime by reducing the number of propagating channels 
in the leads; after the threshold of the lowest transverse mode 
becomes larger than the Fermi energy in the leads, an effec- 
tively ID tunneling barrier is formed in this lowest mode. In 
a MQD, in contrast, tunneling junctions with the leads can be 
made by oxidizing the surface of the dot. The main difference 
for our purposes is that in the case of a MQD there are many 
channels propagating at Ep whereas in SQD's there is only 
one relevant quantum channel. 

By capacitively coupling to a metallic gate, one can control 
the dot's potential, allowing current to flow and bringing out 
mesoscopic effects. 1 For certain values of the gate voltage V g , 
the electrostatic energy difference between N electrons in the 
dot and N+l electrons is balanced by the interaction with the 
gate: an electron can freely jump from the left lead onto the 
dot and then out into the other lead. This process produces a 
peak in the conductance at this V g . By changing the gate volt- 
age one can observe large peaks followed by valleys. Both 
the peak heights and peak spacings fluctuate as one varies V g . 
The distribution and correlations of the peak heights and spac- 
ings have been studied experimentally, as well as theoretically 
using random matrix theory (RMT)illn general there is agree- 



ment between experiment and theory4^ 

Whether the temperature T is larger or smaller than the 
typical width of the resonances T has a large effect on the 
dot's conductance. As long as T ^> T transport proceeds by 
either resonant tunneling (peaks) or the off-resonant process 
known as co-tunneling (valleys). 7 However, when T is suffi- 
ciently small many-body effects become important; in partic- 
ular, at T = if N is odd the conductance can be of order 
the conductance quantum e 2 /h in the Coulomb blockade val- 
leys where naively one expects a strong suppression of the 
conductance. 8 ' 9 Only recently have experiments actually seen 
this enhancementJ^iiii 2 ^ 

The temperature scale at which these many body effects de- 
velop is called the Kondo temperature Tk- 14 This scale is ex- 
ponentially sensitive to the tunneling rates from the leads to 
the dot. Since these tunneling rates depend on the wavefunc- 
tions on the dot, Tk shows mesoscopic fluctuations. 

Experimentally, mesoscopic fluctuations of the Kondo ef- 
fect in SQD's is observed. Perhaps the clearest example is in 
Ref. [l2| where the conductance as a function of magnetic field 
in the Kondo regime is shown. In addition, the frequent obser- 
vation of the Kondo effect in certain valleys but not others 11 
suggests that mesoscopics plays a role. Likewise, the obser- 
vation of the Kondo effect in two adjacent valleys 13 argues for 
the important role of fluctuations^ 

At low temperatures Tr- is the only important energy scale 
in the problem, 14 and thus one can calculate the distribution 
of any property given the distribution of Tk- We will focus 
on the conductance as it is the most relevant experimentally. 

Assuming chaotic dynamics in the dot, we use random ma- 
trix theory to calculate the distribution of Tk in the CB val- 
leys. We restrict our attention to the case of odd N and 
5 = 1/2 and study how the fluctuations depend on the number 
of propagating channels. We go on to calculate the distribu- 
tion of the conductance at T = and discuss the effect of 
finite temperature. 

The Hamiltonian — The Hamiltonian of the system, H = 
Hdot + -ffieads + -Hti consists of the quantum dot's Hamilto- 
nian, the Hamiltonian of the leads, and the tunneling Hamil- 
tonian which describes dot-lead couplings. -ffdot has both a 
one-body part and an interaction part. Impurities and/or scat- 
tering from the boundaries are taken into account statistically 



2 



by using a random matrix model for the one-body Hamilto- 
nian. The most important interactions of the electrons on the 
dot are the charging effect and exchange. 7 Thus, we find 



Hdot = t 
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where Sj is the single-particle energy spectrum on the dot; 
Af is the dimensionless gate voltage used to tune the number 
of electrons on the dot; n is the number operator for excess 
electrons on the dot; S is the spin operator; and Jg and Ec 
are the exchange constant and charging energy, respectively. 

The Hamiltonian of the leads differs in the SQD and MQD 
cases. In the SQD a tunneling junction is made by pinching 
until just after the last propagating mode is cut off - the lead 
is effectively ID. In the MQD, the leads are wide and many 
propagating channels can tunnel through the oxide barrier. In 
the general case of N-^ and Nr channels, we can label all the 
states in the leads by a channel index and ID momentum, 
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where efc + E m are the one-particle energies in the m th chan- 
nel with momentum k. 

Finally, the Hamiltonian for the dot-lead coupling is 
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where t m j are the matrix elements for each of the Ni, + Nr. 
quantum channels tunneling into the j state of the quantum 
dot. \t m j | 2 is proportional to the intensity of the wavefunction 
j in the dot. We assume that it is independent of k in the 
lead since the typical energy scale for changing wavefunctions 
in the clean leads is much larger than other relevant energy 
scales. 

We would like to rewrite this Hamiltonian in terms of en- 
ergy rather than momentum states in the leads. We define new 
operators cj„ ea . = c\ nka \de/ dk^ 1 / 2 that create states with en- 
ergy e = E m + fi 2 k 2 /2m in the m th channel. The normaliza- 
tion is chosen so that [cj,, ct,] = 5(k — k') implies [c e , cL] = 
S(e — e'). In terms of these operators the Hamiltonian is 

Hh=Y^ J dee (CroeaCmetr) 
ma 
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with t m j = t m j \de/ dk\ x / 2 . Note that even for the same e 
the derivatives and the lower limits on the integrals will be 
different for different channels. 

Mapping to single channel Anderson model — The Hamil- 
tonian of the system bears close resemblance to an iV-channel 
impurity problem. However, since we are considering the 
5 = 1/2 odd TV valley Kondo regime, 7 only one of the spa- 
tial states in the dots is of consequence, namely the highest 



energy orbital which we shall refer to as j res . When the am- 
plitude from this one level leaks out through the barriers it is 
intuitively plausible that it couples to only one linear combi- 
nation of the transverse wave-functions in the leads. Hence 
we look for a rotation in the channel basis that chooses the 
correct wavefunction in the leads to couple to level j xes on the 
dot, denoting the new creation and destruction operators in the 
leads z' and z. If we choose 
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where th,R — (X) TO gi, r ^J 1 ' 2 ' on ly one channel couples 
to the quantum dot on each lead. We may choose the remain- 
ing channels in any way just making sure they are orthonormal 
to zi 2- This defines a unitary rotation matrix U the first two 
rows of which we have specified. 15 The Hamiltonian has now 
been effectively reduced to two parts: (1) two single channel 
leads connected to a quantum dot with tunneling matrix ele- 
ments tL and tji, and (2) TVl + TVr — 2 decoupled channels. 16 
The decoupled channels do not contribute, and the solution 
to the first problem is available in the literature: 8 By making 
another rotation in R-L space, one can reduce the problem to a 
1-channel Kondo problem plus a decoupled channel. Thus all 
the thermodynamic properties of the quantum dot are those 
of a one-channel Kondo problem. For the conductance, one 
must rotate back to the original basis of leads. We may use 
the Kubo formula to calculate the conductance at finite tem- 
peratures with the result 17 
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where T a {uS) is the T-matrix of the scattering problem, / is 
the Fermi function, and Go is defined in terms of the level 
widths r m = 2irv m \t mjtea | 2 by 
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So we may write the expression for the conductance as a prod- 
uct of a prefactor Go and a universal function Fk(T /Tk) 
which has been calculated 1819 using the numerical renormal- 
ization group technique. 

Fluctuations of Tk — Because of the exponential sensitiv- 
ity of Tk on the wavefunctions, fluctuations are expected to 
be strong in the weak tunneling regime of a SQD. As one 
increases the number of channels one obtains self-averaging. 
Our goal is to study the crossover between these two limits. 

For a simple Anderson model for the quantum dot, the 
Kondo temperature from a scaling argument 2 ^ is 



T K = (t/r) 1 / 2 exp 
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FIG. 1 : The moments of the distribution of Tk as a function of the 
number of channels in the leads iV c h (Eg. 1121 . The first three mo- 
ments are shown - average (dotted), root-mean-square (dot-dashed), 
and asymmetry (solid). The moments indicate that the distribution 
is badly behaved for small iV c h while it becomes reasonable for 
iVch > 10. 



where eo is negative and measures the energy difference be- 
tween Ep and the singly occupied level on the dot, T is the 
width of the level, and U = 2Ec is the on-site interaction. 
In the case of a quantum dot, we have to correct for the fact 
that the spectrum in the dot is dense, and thus the high energy 
cutoff for the Kondo Hamiltonian is not given by U but by the 
mean spacing A. 7,20 Varying Vg adjusts eo; expressing this as 
a fraction of U by choosing e = — xU(x > 0), we write Tk 
as 



T K = 



A 



1/2 

— j exp 



TTX(1 — x)U 



(10) 



To calculate the distribution of the Kondo temperatures one 
needs an appropriate distribution for T. For simplicity we as- 
sume that the different channels in the lead couple to the dot 
wavefunction in a similar way on average so that the average 
level widths for all the channels are the same although they 
fluctuate independently. We allow, however, for a different 
number of channels in the left and right leads. From random 
matrix theory, it is known that the total level width follows a 
X 2 distribution with N c ^ = (3(Ni,+N-r) degrees of freedom: 4 
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The distribution of tn = Tk\Jx(\ — x)/A follows from 
Eqs. JlOi and ( II Q . While we have not been able to obtain 
a closed-form expression for P(tjf), we can calculate all its 
moments: 
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FIG. 2: Probability density of Tk, P(Tk), calculated numerically 
in four cases: AT c h =1. 2, 4, and 50. The main panels are on a linear 
scale while the insets are the distribution of log Tk- The distributions 
for small 7V c h are very broad while that for N c ^ = 50 is reasonably 
behaved, though not yet Gaussian. Note that applying a magnetic 
field doubles the number of effective channels in the leads - changing 
iV c h = 2 to 4, for instance - and so increases the probability of seeing 
a Kondo effect. We have used T — 0.2(7 and x — 1/2. 



where a = U /2T and K n (z) is the modified Bessel function 
of the second kind. 

Note, first, that when P(T) is broad, as for iV c h = 1 or 2, 
the fluctuations of Tk are on a logarithmic scale and so are 
huge. On the other hand, the Kondo temperature depends on 
only the sum of the decay widths. On adding more channels 
one gets a sharply peaked function for T. Thus, the huge log- 
arithmic fluctuations inherent in the Kondo effect are reduced 
until something well-behaved is obtained. 

In Fig.Q we have plotted the average, rms, as well as cubic 
deviation as a function of the number of channels using our 
analytic formula. The number of channels at which there is a 
crossover, from logarithmic fluctuations to a distribution well 
characterized by its mean and variance, is about /Y c h — 10. 
We also note that the cubic moment and rms decay only alge- 
braically and so the crossover to self-averaging is slow - there 
is no "scale" associated with this crossover. The main features 
of the dependence on 7Y c h is clearly borne out in the numerical 
results for the full distribution shown in Fig. [2] 

Zero Temperature — In the case of many channels coupled 
to the quantum dot, (Tk) is a meaningful quantity and we can 
define a T = regime. To calculate the distribution of con- 
ductance at zero temperature, we use Go = Gk{T = 0) taken 
from Eq. l|8}. This distribution is also a good approximation 
in the case ^> TVr as Tk depends on A^l + and hence 
does not fluctuate as much as Go- 

Using again the random matrix theory result il It for the 
level widths, the above calculation gives a result in closed 
form: 
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FIG. 3: Probability density of the conductance for three different 
sets of (JV L , Nn) in the T = regime. 
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where g is defined by G = 4(2e 2 /h)g. 

We have plotted P(g) in Fig. [3] for three different channel 
realizations to emphasize that the distribution function can be 
changed quite dramatically by changing N-^ and Nr. As ex- 
pected, in the symmetric case the most probable conductance 
is 2(e 2 /h). By putting in some asymmetry we can obtain a 
variety of distributions. Note in particular the highly asym- 
metric case in which the distribution peaks at a conductance 
much smaller than 2(e 2 /h). Since the fluctuations of Tk de- 
pend on the sum of the we expect in a heavily asymmetric 
case that the zero temperature result will show even at finite 
temperatures. Such a situation may be realized in a metallic 



particle when one of the "leads" is an STM tip. 

Conclusions — Our main results are: (1) The Kondo en- 
hanced conductance is given by Q for any number of chan- 
nels in the left and right lead. (2) We have calculated the distri- 
bution function for both Tk (Figs. 1 and 2) and the prefactor 
Go (Fig. 3). (3) At T = 0, fluctuations are dominated by the 
prefactor, and one should look at Jl 31 for the distribution. 

Turning to a comparison with the experiments, we first 
note that most of the experimental dots were not in the deep 
Coulomb blockade regimep2iiiii2ii^ and so application of our 
results is problematic because we have neglected charge fluc- 
tuations. Nonetheless, we use our theory to make a few exper- 
imentally relevant estimates. One of the difficulties in the ex- 
periments is finding several odd valleys in a row which show 
the Kondo effect. 21 The changing characteristics of the dot as 
a function of gate voltages contribute to this difficulty, but the 
mesoscopic fluctuations of Tjc will also play a role. RMT 
predicts that neighboring energy levels have completely un- 
correlated wavefunctions and so that Tk in sequential odd N 
valleys are uncorrected. For the dots of Ref.[T3, for instance, 
we estimate that Tk > T in about 0.4 of the odd valleys (us- 
ing (T/U) = 0.2 and T/A = 0.02). 

Fluctuations of Gk can also be generated by continuous 
tuning parameters such as magnetic field, revealing a corre- 
lation scale. The expected scale is the semiclassical scale 
for changing chaotic wavefunctions. Interestingly, in Ref.fl2l 
these fluctuations are correlated on a scale which is an order of 
magnitude larger. This is consistent with measurements of the 
correlations in Coulomb blockade valleys in the co-tunneling 
regime. 22 Neither measurement is understood at this time. 
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